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Inspired by the obvious discrepancies between experiment and non-relativistic QCD (NRQCD) 
studies of e + e~ — + J/if) + rj c at \/~s ~ 10.6 GeV, we investigate contributions from intermediate 
meson loops as long-range interaction transitions to this process. The intermediate meson loops 
include DD(D*), DD*(D or D*), D*D*(D) and corresponding D 3 intermediate mesons. With the 
constraints from experimental data on the vertex couplings, we find that the intermediate meson 
loops account for 2.7 ~ 6.7 fb of the cross sections within a reasonable range of cut-off energies of 
the factor parameter. We also investigate contributions from the absorptive part and find that it 
accounts for approximately 0.58 ~ 1.38 fb. These results imply that contributions from long-range 
interaction transitions may still play a role in such an energy region. 

PACS numbers: 13.66.Bc, 12.38.Lg, 14.40. Gx 



I. INTRODUCTION 

Recently, one of the most hottest topics is the significant discrepancies between experimental and 
theoretical results for exclusive double-charmonium production e + e~ — > J/ip + r\ c at the center mass 
energy of 10.58 GeV. In 2002, Belle Collaboration firstly reported the exclusive cross section for a[e + e~ — ► 
J/ip + r] c ] x B(j] c ->> 4 charged) = 33±g ±9/6 fl|]. In 2004, Belle Collaboration updated their results 
a[e + e- -> Jjvb + r) c ] x B(r) c ^> 2 charged) = 25.6 ± 2.8 ± 3.4 fb 2] and BABAR Collaboration also 
measured the same quantity, and found a[e + e~ — > J/ip+r] c ] xB(r) c — >> 2 charged) = 17.6±2.8±2.1 fb |3j. 

Theoretically, based on nonrelativistic QCD(NRQCD) factorization approach [4| at leading order (LO) 
in the QCD coupling constant (a s ) and charm-quark relative velocity^) , the predictions given by Braaten 
and Lee Liu, He and Chao [g], and Hagiwara, Kou and Qiao [3] are about 2.3 ~ 5.3 fb, which are 
about one order of magnitude smaller than the experimental measurements. In order to solve this puzzle, 
many solutions have been proposed. In Ref. [1] , the authors consider the corrections of next-to-leading 
order (NLO) in a s , which enhanced the cross section with a K factor (the ratio of LO plus NLO to LO) 
of about 1.8 ~ 2.1. Authors of Ref. [5J find that including relativistic corrections (order of v 2 ) will lead to 
the K ~ 2.0l 1 ° 1 9 , but with large uncertainties already bared with the NRQCD matrix elements. In order 
to reduce the uncertainties of the order-u 2 matrix elements, a potential model is applied to calculate the 
quarkonium wave function [9( , by which the relativistic correlations and corrections of NLO in a s can be 
treated [l(|, and a cross section of 17.5 ±5.7 fb is obtained. In Ref. [ll|, the authors present their results 
by resummation of relativistic corrections which contains several refinements, such as nonperturbative 
NRQCD matrix elements, inclusion of the effects of the running of a s etc. They conclude that the 
discrepancies between the theoretical and experimental can be understood. By determining the matrix 
elements in a different way, Ref. pj| obtains a value of 20.04 fb. A complete computation of leading and 
NLO contributions are recently presented by Gong and Wang [13| |. which is a directly confirmation of 
results from Ref. [8[ . Treatments using light-front approach such as [l4| [HI, [l6|, [TtJ also provide other 
possible solutions. 

In brief, the present theoretical studies show that NLO contributions turn out to be important in this 
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double-charmonium production, which simultaneously raises concerns about the perturbation expansion. 
On the other hand, large contributions from NLO corrections suggest that nonperturbative mechanisms 
may start to play a role. This corresponds to the long-range part of the strong interactions, among 
which intermediate meson loops could be a natural explanation for the cross section enhancement in 
e+e- -> J/ip + f] c at ,/i ~ 10.6 GeV. 

Intermediate meson loop (IML), or intermediate meson rescattering (or sometimes it is presented as 
final state interaction (FSI) with on-shell approximation), as one of the important non-perturbative 
transition has been investigated extensively in heavy meson decays. Cheng et al. studied the long- 
distance rescattering effects in hadronic B decays [18]. From the data accumulated at B factories and 
CLEO, it was found that soft final-state rescattering effects played an essential role in B physics. In 
Ref. [20| | . Liu et al. also found that contributions from the hadronic loops turned to be important in 
charmonium hadronic decays. Recently, a systematic investigation of the IML effects in quarkonium 
hadronic decays involving OZI-rule violations and isospin breaking reveals that IML plays a role in many 
places and sometimes can even compete against direct scatterings [U [2^, HH . 

In e + e~ — > J/ip + r) c , the absorptive contributions via IML are supposed to be present, and they do not 
interfere with the LO NRQCD transitions. With constraints from experiment on the effective couplings, 
we can explore the IML mechanism in e + e~ — ► J /if) + rj c . 

The paper is organized as follows: In section II, we describe the IML model with effective Lagrangians. 
The numerical results are given in Section III. Discussion and summary are given in Section IV. 

II. THE MODEL 
A. Effective Lagrangians 

In reaction e + e~ — > J /ip + rj c , the final J/ip and r\ c can be produced by direct production of two cc 
pairs. This process is taken care of by pQCD transitions. In our model, we consider indirect production 
process where J/ijj and rj c are produced by the intermediate meson loops as shown by Fig. [TJ 

In principle, we should include all the possible intermediate meson loops in the calculation. But, 
the break-down of the local quark-hadron duality allows us to pick up the leading contributions as a 
reasonable approximation in the practical calculation [H], Hfl . Some possible leading IML contributions 
are shown in Fig. [2] Since those intermediate mesons can be off-shell at s/s ~ 10.6 GeV, apart from the 
absorptive part of the loop transition amplitudes, the real part can also contribute. This is different from 
the FSI approach where only the absorptive amplitudes are considered. In this work, we shall investigate 
both. 

In order to evaluate the diagrams, we adopt the following effective Lagrangians: 
£jdd = g^DoiDd^D - d ll DD}A^ 
£ lD r>* = -ig lD D^ a p^d a Al 3 d^D* v D + h.c., 
t-tD-D* = g jD . B .{A»(d tl D"'D* v -D*''d M D* v ) 

+ {d^A v D* u - A^D* V )D^ + D*»(A»d^D\ - d^D*")}, 

^dd« = -ig^ D D*e a ^ v d a ^d»D* u D + h.c. (1) 

where A^ 1 and tp^ are photon and J/ip vector-meson fields, respectively, and D is the pseudoscalar- 
meson field with D = (D°,D+,Df) and D = (D° , D~ , D~) T '; D* is the vector-meson field with D* = 
(D*°, D* + , D* s + ) and D = [D* , D*~ , D* S ~) T . The coupling constants appearing in the above equations 
will be determined as follows. 

For the couplings of photon with D-meson, such as g^oDi i 9~{D*d* i we can extracted them from 

the experimental data [26j. Applying the effective Lagrangians, the 7 couplings can be obtained: 

2a e \Pd d 
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9-/DD* 



16s 5 / 2 cr e+e -_ Dj5 » p /2 

a e \P D >\ C DD* 



_ f ^ s5 ^ 2cr e+e--^D'D* m D^ l/2 ( r )\ 

9^D>D* - I B , 1 , W 



with 



C D d> = ™d* - 2 ( m b + s)m%, + {m 2 D - s f + -s\P 3 . | 2 

C D , D , = s 3 + 6s 2 m 2 D , -28sm D . - 48m 6 D , - -{12m%, - 6sm 2 D , +s 2 )|P^.| 2 , (3) 

where Pjj, Pd», and P' D , are the final state three-momentum for the processes e + e~ — > DD, DD*, 
D*D* , respectively; a e = 1/137 is the fine structure-constant; a are the corresponding cross sections. In 
the SU(3) limit, we have the following relation g lDs u 3 = 9 iD d, 5 7 d s d* = 9 7 dd« and 5 7 d s *d s * = 9~ ( d*d*- 
The values are listes in Table HI 

For the J/ ipDD* coupling, we use the relation in the heavy quark mass limit [U HI] , 

9i,DD' = 9iPDd/Md, (4) 

where Mo is the mass scale of the D/D* mesons. We adopt the coupling constants g^uB = 7-44 and 
9<I>dd* = 3.84 GeV" 1 from Ref. [2^ . and g Va D*D* = 9ibDDi 9ij a DD* — 9^dB* m the following calculation. 



B. Intermediate meson loops contribution 

We take DD(D*), DD*(D or D*), D*D*(D) and corresponding D s intermediate mesons into account. 
The transition amplitude for e + e~ — > J/%p + rj c via intermediate meson loops can be expressed as follows: 

Mfr = v^m-ie^u^e)^ I P^^^HpII (5) 
1 s J (2tt) 4 axa 2 a 3 

where Ti^,3 are the vertex functions of which the detailed expressions have been given in Refs. [HI, HI]. 
We include them in the Appendix for all the above mentioned loops. Four-vector momentum p 2 is for 
the exchange meson. Variables ai = p\ — to 2 , a 2 = p 2 — m 2 ,, and 03 = p\ — m 2 are the denominators of 
the propagators of the intermediate mesons, respectively. 

In the above equation a form factor T(p 2 ) is introduced to take into account the off-shell effects of the 
exchanged mesons and also kill the divergence of the integrals. The following commonly-used form factor 
is adopted: 

(«> 

where n — 0, 1, 2 correspond to different treatments of the loop integrals. A is the cutoff energy, which 
should not be far away form the physical mass of the exchanged particles. Since there are different 
particles exchanged in the meson loops, a useful parameterization is as follows (l8| : 

A = m ex + uAqcd, (7) 

where Aqcd = 220 MeV and a is a tunable parameter; m ex is the mass of the exchanged meson. This 
way of parameterizing the cut-off energies is slightly different from the previous works [U [22|, [23| , where 
A was universal for all the meson loops. Here, the mass differences between the exchanged particles are 
taken into account. The parameter a will then be constrained by experimental data. In this work, we 
apply the dipole form factor in the calculation, i.e. n = 2. The relation between the loop amplitudes 
with the dipole and monopole form factors are given in Appendix III. 
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C. Absorptive contributions from the IML 

It is important to examine the absorptive part of the IML transitions. Non-negligible contributions 
from the absorptive part would be a strong evidence for possible long-range transitions in the double- 
charmonium productions. 

Taking D(pi) (or D*(pi)) and D(p 3 ) (or D*{p%)) as on-shell particles, we obtain the absorptive tran- 
sition: 

1 ie p 



2 yre " ,H ' s 

d 3 pi d 3 p 3 4 4 , TiT 2 T 3 ., 

fO ^30F In ^30F ( 2?F 6 (Pe+Pe-Pl-P3)— 2~T~- T^fe) 

(2Tty2k,\ {2-K) 6 2hj3 p^ — m,2 + irriiV 2 



r 1 |pi|dcos0 r!r 2 T 3 

s 7_! 16^^ — m| + «m 2 r 2 

where is the angle between p\ and p^. The kinematic definitions have been given in Appendix II. The 
vertex functions 21,2,3 and form factor T(j$£) are the same as the previous definitions. 

D. EM transition via T(4S) intermediate meson 

The present data from Belle were taken at y/s = 10.58 GeV and 10.52 GeV, respectively. The higher 
one corresponds to the mass of T(45) while the lower one is a side-band measurement. However, the 
datum samples are not sufficient for determining the cross section differences between these two energies. 
Therefore, the role played by T(4S') from the present Belle results [H, 0, HfjJ is still unclear, and the 
T(45) — > J/i>+r) c coupling is also unknown. In Ref. [29|, exclusive decays of T(45) to double-charmonium 
states were evaluated and the branching ratio for T(45) — > J/tp + rj c was predicted to be at an order of 
10 -9 . At amplitude level, such a small contribution can still produce some structures at the T(45) mass 
due to the destructive interference between the resonance strong decay and continuum amplitudes. 

In our calculation we take into account the contribution of the T(45). Its interference with the IML 
amplitudes is useful for us to examine the model-dependence of the IML transition and sensitivities of 
the cross sections to its coupling to J/ipr] c . 

Fortunately, the EM coupling for T(45) — * e + e~ has been relatively well measured |3Cb| . We can then 
determine the V^* coupling eMy/fv by the vector meson dominance (VMD) model |3l| . 



e 



3r, 



2a P 



1/2 



(9) 



where Ty_, e + e - is the partial decay width, |p e | is the electron three-momentum in the vector meson rest 
frame, and a e = 1/137 is the fine-structure constant. It should be noted that this form of interaction 
is only an approximation and can have large off-shell effects arising from either off-shell vector meson or 
virtual photon fields. In this approach we consider such effects in the VjP coupling form factor which 
will then be absorbed into the energy-dependent widths of the vector mesons. 
The transition amplitude corresponds to Fig. [3] is 

M fi = «^ VWKWj^^f ■ 

In case that the partial width of rx^^ c is available, we can extract the Tipr) c coupling via 

ffT^JPy-l"' , 1(V . 

1 T->V"7c - 77, > l iU J 

1Z7T 



where |p^,| is the three-momentum of J/ijj in the T meson rest frame. 
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The total transition amplitudes apart from the leading QCD contribution can thus be expressed as 

M fi = Mf t M + Mf° op . (11) 
We estimate the differential cross section via: 

B-akigfi5>*r. (12) 

spin 

where p e is the three-momentum of the initial electron (positron) in the overall cm. system. The mass 
of the electron has been neglected, i.e. p e = E cm /2 is applied. 



III. NUMERICAL RESULTS 



As mentioned before, the present experimental data at y/s = 10.52 and 10.58 GeV cannot tell the 
role played by T(4S f ). But we can still define quantity R, which is the e + e~ — > J/%p + r\ c cross section 
ratio at y/s = 10.52 GeV to that at y/s — 10.58 GeV, i.e. R = ctio^/ctio.ss- It is a function of both 
resonance direct transition and IML amplitudes, and we expect that it should not be sensitive to the cut- 
off energy introduced to the IML though the exclusive IML contributions may have strong dependence 
on it. Unfortunately, no information about the T(45 l ) coupling to J/ipr] c is available. We then simply 
assume R = 0.5 at a = 1.8 to fix the coupling gr4>n c = 5.26 x 10~ 5 GeV -1 . This corresponds to 
BR(T — > J/tprjc) = 2.9 x 10~ 7 . One should not take this value seriously since our strategy here is to fix 
gxipr/c i and then examine the sensitivity of R to the form factor parameter a as a test of the behavior of 
the IML. 

In Fig. HI by fixing g~t^n c = 5.26 x 10 -5 GeV -1 , we plot the ratio R with a varying a. It shows that 
within the commonly accepted range of a = 1.6 ~ 2.0, the ratio R varies from 0.43 to 0.57. This is an 
indication of insensitivity of the loop contributions to the form factors within a reasonable range of the 
cut-off energy. We also find that such a property is retained even for larger gr^, Vc . In such a case, the cross 
section at the mass of Y(4S) is enhanced by the resonance contribution while the side-band cross section 
is relatively small. A quantitative determination of the T(4S") requires future precise measurements of 
the cross sections at both y/s — 10.52 and 10.58 GeV. 

We also plot the cross-section-dependence on a at y/s — 10.52 and 10.58 GeV in Fig. [SJ It shows that 
these two cross sections increase with a slowly and also appear to be stable. Again the difference between 
these two cross sections are due to the assumed contribution from T(45 f ). 

In Fig. O we fix the form factor parameter a = 1.8 and then investigate the effects from the T(45) 
at different R values. The total cross sections from the IML plus resonance T(45) are presented. A 
smaller value of R corresponds to a larger gTipr tc coupling in the present choice of a constructive relative 
phase. It is possible that their interference is destructive and a dip will appear at y/s — 10.58 GeV [29j]. 
In this sense, we emphasize again that it is essential to have precise data for the cross sections at both 
y/s = 10.52 and 10.58 GeV. 

In Fig. [71 we plot the total cross section again, but with several different a values. The solid curve is 
the same as that in Fig. [6] with a — 1.8. Again, some sensitivities to a are highlighted. We also present 
the results without T(45), i.e. exclusive cross sections from the IML, as denoted by the dot-dashed curve. 
A flat behavior is observed along with y/s. 

It is interesting to see that the total cross section without the T(45) has the same order of magnitude 
as the pQCD leading order results. In Tab. HTlwe list the total cross sections with different a values. In a 
range of a = 1.6 ~ 2.0 the total cross sections increase from 2.74 fb to 6.76 fb. Eventually, this range is 
acceptable for the form factor uncertainties which is unavoidable in this effective Lagrangian approach. 
Similarly, the side-band cross sections at y/s = 10.52 GeV (i.e. can be viewed as no T(45) contributions) 
vary from 1.3 fb to 3.5 fb. 

In order to further clarify the role played by the IML, we calculate the absorptive contributions from 
the loops by making the on-shell approximation for the intermediate mesons. In Fig. El the a-dependence 
of the cross section at y/s = 10.58 GeV is presented. Since the absorptive part will be present in the 
transition amplitudes, its stability within the commonly accepted range of the form factor parameter 
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turns to be essential for understanding its role in this process. As shown by Fig. [HI with a from 1.6 - 2.0, 
the absorptive cross section varies from ~ 0.6 — 1.4 fb, and can be regarded as quite stable. We also list 
the cross sections with different values of a in Tab.[H]to compare with the full loop calculations. 

We then examine the total absorptive cross sections with fixed a values in Fig. [9l It shows that with 
a = 1.8, the cross section is about 0.95 fb, and its energy-dependence is rather weak. Again, larger values 
for a produce larger cross sections. The most interesting feature is that the absorptive cross sections are 
nearly the same order of magnitude as the LO results of NRQCD. 

One should be cautioned in the understanding of the result of Fig. [9l It eventually reflects the energy- 
dependence of the IML form factors, and with the couplings of e + e~ — * DD etc fixed at y/s ~ 10.56 GeV. 
In this approach QCD hard interactions have been contained in those couplings derived at y/s ~ 10.56 
GeV, i.e. g lDD , g 1 DD' an d 5 7 _d*d*- The form factor in Eq. ([5|) will then take care of the soft interactions 
arising from the intermediate meson loops. Hence, one should not compare the energy-dependence of 
Fig. [5] with the QCD factorization (g, 33] for the direct production of double charmonia. 

Also, it should be pointed that the IML contributions as a source of long-range interaction transition are 
different from the soft exchanges in the pQCD factorization 33] . It is not a full correspondence of the soft 
QCD. Ideally, we expect that the sum of all the possible intermediate meson loops be equivalent to the soft 
exchanges in the pQCD factorization according to the quark-hadron duality argument [U [H[ . In reality 
the break-down of the local quark-hadron duality allows us to pick up the leading IML contributions as 
a first-order approximation. 

IV. SUMMARY 

We estimate the intermediate meson loop contributions to e + e~ — > J/ip + r) c in an effective Lagrangian 
theory. By applying the available experimental information to the constraints of the meson-meson cou- 
pling vertices, we investigate contributions from DD{D*), DD*(D), DD*{D*) and D*D*(D). The 
model-dependence mainly comes from the form factors adopted for the loop integrals. Fortunately, it 
shows that the results do not vary dramatically with the cut-off energies within the commonly accepted 
range. The resonance T(45) effects are also estimated. Precise measurements of the cross sections at the 
resonance energy and side-band will be able to provide useful information on the T(45). 

In order to clarify the IML contributions, we also calculate the absorptive part of the loops and find 
that the IML contributes nearly the same order of magnitude as the LO of NRQCD. It is likely that the 
long-range IML be an important mechanism apart from the NRQCD LO transitions. Qualitatively, the 
long-range IML transitions may have some overlaps with the NLO processes if they are not obviously 
suppressed. Note that the IML contributions go to the effective coupling of the anti-symmetric tensor 
at hadronic level. It suggests that the relativistic corrections may also have some overlaps with the IML 
mechanism. To understand the large cross sections from the experiment, one perhaps should consider 
both short and long-range transitions to obtain an overall consistent prescription. Again, we emphasize 
the importance of extracting the precise cross sections at and off resonance T(4S f ). 
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Appendix I: Vertex functions for the intermediate meson loops 

For DD(D*), the vertex functions are 

T\ = igi(pi - pa) ■ e 7 

T 2 = i92£ a flvvV%£%V2 e 2 ( 13 ) 

T 3 = i93(Pr,a +P3) ■ £ 2 

where gi, and g 3 are the coupling constants at the meson interaction vertices (see Fig. da)). The four 
vectors, Pjm, and p Vc are the momenta for the final state J ftp and r\ c meson; The four- vector momentum, 
Pi, P2, and P3 are the intermediate mesons, respectively. 

As shown by Fig. |TJ the vertex functions for the DD*(D) + c.c. loop are 

T 2 = i/ 2 (P2 - Pi) ■ H > ( 14 ) 
T 3 = ihij>vc -P2) ■ £3 • 

where /i,2,3 are the coupling constants. 

We also consider the transition amplitude from the intermediate DD*{D*) + c.c. loop , which can be 
expressed as the following formula: 

Ti = ihxEap^p^sPp^e^ , 

T 2 = ih 2 £ a '0'n'v'P2 e 2 Pi, £ i > ( 15 ) 
T 3 = i/l 3 £a"/3" A1 'V"P2" £ 2 P3 £ 3," 

where /ii,2,3 are the coupling constants. 

Transition amplitude from the intermediate D* D*(D) + c.c. loop can be writen as Eq. (|16p 

Ti = iAi [-e 7 • (pi - p 3 )ei • e 3 + 2pi • s 3 si ■ e 7 + 2ei • p 3 e 1 ■ £3] 
T 2 ee iXiSap^pfefp^e^ (16) 
T 3 = iA 3 (p,, c -p 2 ) -£3 

where Ai^.3 are the coupling constants. 



Appendix II: Four-vector momentum expressions and wave functions 



In order to calculate the absorptive amplitude of e + e~ — > J /ip + rj c , it is convenient if we choose 
the overall cm. frame and choose the z axis along the three-vector momentum of J/ip. After taking 
m e = m e + ~ 0, we have 



and 



P% 
Pi 

p\ 

Pi 

a 

Ps 

Pt 

I A* 

Pe 



(^,0,0, \p f \) T 

= (Ei , I pi I sin 9 cos <p, \pi | sin 9 sin ip, \pi | cos ( 

= (E 2 ,- \pi I sin 9 cos <p, - |pi | sin 6> sin </?, \p f \ ~ \p x 
= (E 3 ,—\pi \ sm9 cosip,—\pi \ sin 9 sin <p, —\pi\ cos 9) 
= (E,E sin 9\ cos ipi , £7 sin 9\ sin 931, Ecos9\) T 

= (E, —Esin9i cost/?i, —Esm9i sint^i, —Ecos9i) 7 



) T 

cost 



<(±1) 

<(0) 

U S (Pe) 



T^(0,l,±z,0) T 
(|^ / |/M^ > O J O ) B /M^) T 
77|(x( S ) T ,0,0r 
=4(0,0,x(*) T ) T 



(17) 



(18) 
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where x( s ) is the wave function of the electron and positron; {9\,ip\) are the azimuth angles opened by 
the three-vector momentum. 

1 / COS % \ , 1, /-sin%e-^\ f . 



Appendix III: Useful formula 



Started with the case of no form factor in the integral, Eq. ([5|) can be expressed as the power of the 
four- vector momentum of the exchange- meson (p 2 ), and it is easy to shown that the power of p 2 is no 
more than three. We apply Feynman parameter dimensional regularization scheme to express the integral 
as a linear combination of the following forms: 



d^ V 2 










f dxdy- 


(27T) 4 [(^ 


"Pip) 2 - 


m 2 ][p 2 -TO 2 ][(p 2 - 




16tt 2 J 


d 4 p 2 










1 dxdy 


(27T) 4 [(pa 


-P^Y - 


m 2 ][p 2 - m 2 2 ][{p2 - 


^P 2 ) ~ m 2 } 


16tt 2 1 


d 4 P 2 




V%vIpI 






1 dxdy 


(27T) 4 [(P2 


-Pv) 2 ~ 


m 2 ][p 2 - m 2 ][{p 2 - 


|-p2) _ TO 2] 


16tt 2 L 



M 2 -P 2 ' 



m 2 -p 2 4 V J 

pi± pv per 



M 2 - P 2 



- -T&iS^P" + 6 titT P" +8 v ' T P fi )], (20) 

where P = yp v — xp^p, M 2 = x(p 2 b — m 2 ) + y(p 2 — m 2 ) — (1 — x — y)m 2 ; e is a infinitesimal parameter. 
The details can be found in Refs. [34l.l35l|. 

To cancel the divergence with a monopole form factor, we just replace the exchanged-meson mass by 
the cut-off energy A to obtain the transition amplitude since the divergent term is independent of the 
intermediate meson masses: 

MJ°{m l7 m2,m 3 ,A) = Mfi{m\, ra 2 , m 3 ) - Mfi(mi, A, ra 3 ), 

where M. /i(mi, m 2 , m 3 ) denotes the amplitude of Eq. ^ with T{p 2 ) = 1. 
For the case of a dipole form factor, we use the identity 

A 2 -m 2 \ 2 , A 2 — m 2 (A + 5) 2 — m 2 , . 

v _ 1 -_ y ,- / ~M A2 „2 /A , ^2 » ( 21 ) 

to have 



A 2 - p 2 J MA 2 -p 2 (A + <S) 2 -p 2 ' 



A^7 4 (mi,m 2 ,m 3 , A) = M^°{m 1: m 2 , m 3 , A) + lim f— M^ (m 1; A, m 3 , A + 5) 

5 — s-0 zAo 
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coupling constants 


value 


Cross section (pb) 


9-yDD 


4.81 xlO" 3 


0.04 




2.73 xW' 6 GeV _i 


0.71 




1.10 xlO" 3 


0.65 



TABLE I: The coupling constants of photon with D-meson determined in e e 
sections are taken from Ref. 12611. 



DD, DD", D*D* . The cross 



a 


1.6 


1.7 


1.8 


1.9 


2.0 


O"10.58(/k) 


2.74 


3.52 


4.43 


5.51 


6.76 


"10.58 (/&) 


1.28 


1.64 


2.06 


2.55 


3.11 




0.587 


0.741 


0.924 


1.14 


1.38 



TABLE II: e + e" — > J /ip + rj c cross sections with different a values at y^s = 10.58 GeV. The ratio R — 0.5 is 
applied, a is the results from the IML with T(4S) contribution^' without T(4S) contribution), while a a is from 
the absorptive transitions where no T(4S) considered. 
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FIG. 3: The J/lp + rj c production via resonance T(45). 
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FIG. 4: The sensitivity of R to the form factor parameter a with the TJiprj c coupling fixed at R = 0.5 and 
a = 1.8. 
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FIG. 5: Sensitivities of the e + e — > J/ip + r\ c cross section to the form factor parameter a at yfs = 10.58 (solid 
line) and y/a = 10.52 GeV (dashed line) with R = 0.5. 
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FIG. 6: The y/s evolution of e + e — > J/ip + Tj c cross section with different R values which correspond to different 
branch ratios for T(4s) — > ip + r\ c . The form factor parameter a — 1.8 is adopted. 




FIG. 7: The y/s evolution of e + e — ■> J/ip + rj c cross section with R = 0.5 but varying values of a = 1.7 (dashed 
line), 1.8 (solid line) and 1.9 (dotted line). The dot-dashed line denotes the result without T(4S') contribution at 
a = 1.8. 
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FIG. 8: The Q-dependence of the absorptive cross sections from the meson loops. 
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FIG. 9: The energy evolution of the absorptive cross sections from the meson loops with different a values. 



